Search for r\' decays into vector meson pairs 



M. Ablikim 1 , M. N. Achasov 5 , D. Alberto 40 , F. F. An 1 , Q. An 38 , Z. H. An 1 , J. Z. Bai 1 , R. Baldini 19 , Y. Ban 25 , 
J. Becker 2 , N. Berger 1 , M. Bertani 19 , J. M. Bian 1 , E. Boger 17a , 0. Bondarenko 18 , I. Boyko 17 , R. A. Briere 3 , 
V. Bytev 17 , X. Cai 1 , A. C. Calcaterra 19 , G. F. Cao 1 , J. F. Chang 1 , G. Chelkov 17a , G. Chen 1 , H. S. Chen 1 , 
J. C. Chen 1 , M. L. Chen 1 , S. J. Chen 23 , Y. Chen 1 , Y. B. Chen 1 , H. P. Cheng 13 , Y. P. Chu 1 , D. Cronin-Hennessy 37 , 
H. L. Dai 1 , J. P. Dai 1 , D. Dedovich 17 , Z. Y. Deng 1 , I. Denysenko 176 , M. Destefanis 40 , Y. Ding 21 , L. Y. Dong 1 , 
M. Y. Dong 1 , S. X. Du 43 , J. Fang 1 , S. S. Fang 1 , C. Q. Feng 38 , C. D. Fu 1 , J. L. Fu 23 , Y. Gao 34 , C. Ceng 38 , 
K. Goetzen 7 , W. X. Gong 1 , M. Greco 40 , M. H. Gu 1 , Y. T. Gu 9 , Y. H. Guan 6 , A. Q. Guo 24 , L. B. Guo 22 , Y.P Guo 24 , 

Y. L. Han 1 , X. Q. Hao 1 , F. A. Harris 36 , K. L. He 1 , M. He 1 , Z. Y. He 24 , Y. K. Hcng 1 , Z. L. Hou 1 , H. M. Hu 1 , 
J. F. Hu 6 , T. Hu 1 , B. Huang 1 , G. M. Huang 14 , J. S. Huang 11 , X. T. Huang 27 , Y. P. Huang 1 , T. Hussain 39 , C. S. Ji 38 , 
Q. Ji 1 , X. B. Ji 1 , X. L. Ji 1 , L. K. Jia 1 , L. L. Jiang 1 , X. S. Jiang 1 , J. B. Jiao 27 , Z. Jiao 13 , D. P. Jin 1 , S. Jin 1 , 
F. F. Jing 34 , N. Kalantar-Nayestanaki 18 , M. Kavatsyuk 18 , W. Kuehn 35 , W. Lai 1 , J. S. Lange 35 , J. K. C. Leung 33 , 
C. H. Li 1 , Cheng Li 38 , Cui Li 38 , D. M. Li 43 , F. Li 1 , G. Li 1 , H. B. Li 1 , J. C. Li 1 , K. Li 10 , Lei Li 1 , N. B. Li 22 , 
Q. J. Li 1 , S. L. Li 1 , W. D. Li 1 , W. G. Li 1 , X. L. Li 27 , X. N. Li 1 , X. Q. Li 24 , X. R. Li 26 , Z. B. Li 31 , H. Liang 38 , 
Y. F. Liang 29 , Y. T. Liang 35 , X. T. Liao 1 , B. J. Liu 32 , C. L. Liu 3 , C. X. Liu 1 , C. Y. Liu 1 , F. H. Liu 28 , Fang Liu 1 , 
Feng Liu 14 , H. Liu 1 , H. B. Liu 6 , H. H. Liu 12 , H. M. Liu 1 , H. W. Liu 1 , J. P. Liu 41 , K. Liu 25 , K. Liu 6 , K. Y. Liu 21 , 
Q. Liu 36 , S. B. Liu 38 , X. Liu 20 , X. H. Liu 1 , Y. B. Liu 24 , Y. W. Liu 38 , Yong Liu 1 , Z. A. Liu 1 , Zhiqiang Liu 1 , 
Zhiqing Liu 1 , H. Loehner 18 , G. R. Lu 11 , H. J. Lu 13 , J. G. Lu 1 , Q. W. Lu 28 , X. R. Lu 6 , Y. P. Lu 1 , C. L. Luo 22 , 
M. X. Luo 42 , T. Luo 36 , X. L. Luo 1 , M. Lv 1 , C. L. Ma 6 , F. C. Ma 21 , H. L. Ma 1 , Q. M. Ma 1 , S. Ma 1 , T. Ma 1 , X. Ma 1 , 
X. Y. Ma 1 , M. Maggiora 40 , Q. A. Malik 39 , H. Mao 1 , Y. J. Mao 25 , Z. P. Mao 1 , J. G. Messchendorp 18 , J. Min 1 , 
T. J. Min 1 , R. E. Mitchell 16 , X. H. Mo 1 , N. Yu. Muchnoi 5 , Y. Nefedov 17 , I. B. Nikolaev 5 , Z. Ning 1 , S. L. Olsen 26 , 

Q. Ouyang 1 , S. Pacetti 19 , J. W. Park 26 , M. Pelizaeus 36 , K. Peters 7 , J. L. Ping 22 , R. G. Ping 1 , R. Poling 37 , 
C. S. J. Pun 33 , M. Qi 23 , S. Qian 1 , C. F. Qiao 6 , X. S. Qin 1 , J. F. Qiu 1 , K. H. Rashid 39 , G. Rong 1 , X. D. Ruan 9 , 
A. Sarantsev 17c , J. Schulze 2 , M. Shao 38 , C. P. Shen 36d , X. Y. Shen 1 , H. Y. Sheng 1 , M. R. Shepherd 16 , X. Y. Song 1 , 
S. Spataro 40 , B. Spruck 35 , D. H. Sun 1 , G. X. Sun 1 , J. F. Sun 11 , S. S. Sun 1 , X. D. Sun 1 , Y. J. Sun 38 , Y. Z. Sun 1 , 
Z. J. Sun 1 , Z. T. Sun 38 , C. J. Tang 29 , X. Tang 1 , H. L. Tian 1 , D. Toth 37 , G. S. Varner 36 , B. Wang 9 , B. Q. Wang 25 , 
K. Wang 1 , L. L. Wang 4 , L. S. Wang 1 , M. Wang 27 , P. Wang 1 , P. L. Wang 1 , Q. Wang 1 , Q. J. Wang 1 , S. G. Wang 25 , 

X. L. Wang 38 , Y. D. Wang 38 , Y. F. Wang 1 , Y. Q. Wang 27 , Z. Wang 1 , Z. G. Wang 1 , Z. Y. Wang 1 , D. H. Wei 8 , 
Q.G. Wen 38 , S. P. Wen 1 , U. Wiedner 2 , L. H. Wu 1 , N. Wu\ W. Wu 21 , Z. Wu\ Z. J. Xiao 22 , Y. G. Xie\ Q. L. Xiu 1 , 

G. F. Xu 1 , G. M. Xu 25 , H. Xu 1 , Q. J. Xu 10 , X. P. Xu 30 , Y. Xu 24 , Z. R. Xu 38 , Z. Z. Xu 38 , Z. Xue 1 , L. Yan 38 , 
W. B. Yan 38 , Y. H. Yan 15 , H. X. Yang 1 , T. Yang 9 , Y. Yang 14 , Y. X. Yang 8 , H. Ye 1 , M. Ye 1 , M.H. Ye 4 , B. X. Yu 1 , 
C. X. Yu 24 , S. P. Yu 27 , C. Z. Yuan 1 , W. L. Yuan 22 , Y. Yuan 1 , A. A. Zafar 39 , A. Zallo 19 , Y. Zeng 15 , B. X. Zhang 1 , 
B. Y. Zhang 1 , C. Zhang 23 , C. C. Zhang 1 , D. H. Zhang 1 , H. H. Zhang 31 , H. Y. Zhang 1 , J. Zhang 22 , J. Q. Zhang 1 , 
J. W. Zhang 1 , J. Y. Zhang 1 , J. Z. Zhang 1 , L. Zhang 23 , S. H. Zhang 1 , T. R. Zhang 22 , X. J. Zhang 1 , X. Y. Zhang 27 , 
Y. Zhang 1 , Y. H. Zhang 1 , Y. S. Zhang 9 , Z. P. Zhang 38 , Z. Y. Zhang 41 , G. Zhao 1 , H. S. Zhao 1 , Jiawei Zhao 38 , 
Jingwei Zhao 1 , Lei Zhao 38 , Ling Zhao 1 , M. G. Zhao 24 , Q. Zhao 1 , S. J. Zhao 43 , T. C. Zhao 1 , X. H. Zhao 23 , 
Y. B. Zhao 1 , Z. G. Zhao 38 , Z. L. Zhao 9 , A. Zhemchugov 17a , B. Zheng 1 , J. P. Zheng 1 , Y. H. Zheng 6 , Z. P. Zheng 1 , 
B. Zhong 1 , J. Zhong 2 , L. Zhong 34 , L. Zhou 1 , X. K. Zhou 6 , X. R. Zhou 38 , C. Zhu 1 , K. Zhu 1 , K. J. Zhu 1 , 
S. H. Zhu 1 , X. L. Zhu 34 , X. W. Zhu 1 , Y. S. Zhu 1 , Z. A. Zhu 1 , J. Zhuang 1 , B. S. Zou 1 , J. H. Zou 1 , J. X. Zuo 1 

(BESIII Collaboration) 

1 Institute of High Energy Physics, Beijing 100049, P. R. China 
2 Bochum Ruhr-University, 44^80 Bochum, Germany 
3 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
4 China Center of Advanced Science and Technology, Beijing 100190, P. R. China 
5 G.I. Budker Institute of Nuclear Physics SB RAS (BINP), Novosibirsk 630090, Russia 

6 Graduate University of Chinese Academy of Sciences, Beijing 100049, P. R. China 

7 GSI Helmholtzcentre for Heavy Ion Research GmbH, D-64291 Darmstadt, Germany 

8 Guangxi Normal University, Guilin 541004, P. R- China 
9 Guangxi University, Naning 530004, P- R- China 
10 Hangzhou Normal University, Hangzhou 310036, P. R. China 
11 Henan Normal University, Xinxiang 453007, P. R. China 
12 Henan University of Science and Technology, Luoyang 471003, P. R. China 



2 



13 Huangshan College, Huangshan 245000, P. R. China 
14 Huazhong Normal University, Wuhan 430079, P. R. China 
15 Hunan University, Changsha 410082, P. R. China 
16 Indiana University, Bloomington, Indiana 47405, USA 
17 Joint Institute for Nuclear Research, 141980 Dubna, Russia 
18 K VI/ University of Groningen, 9747 AA Groningen, The Netherlands 
19 Laboratori Nazionali di Frascati - INFN, 00044 Frascati, Italy 
20 Lanzhou University, Lanzhou 730000, P. R. China 
21 Liaoning University, Shenyang 110036, P. R. China 
22 Nanjing Normal University, Nanjing 210046, P. R. China 

23 Nanjing University, Nanjing 210093, P. R. China 

24 Nankai University, Tianjin 300071, P. R. China 

25 Peking University, Beijing 100871, P. R. China 

26 Seoul National University, Seoul, 151-747 Korea 

27 Shandong University, Jinan 250100, P. R. China 
Shanxi University, Taiyuan 030006, P. R. China 

Sichuan University, Chengdu 610064, P. R- China 
30 Soochow University, Suzhou 215006, P. R. China 
31 Sun Yat-Sen University, Guangzhou 510275, P. R. China 
32 The Chinese University of Hong Kong, Shatin, N.T., Hong Kong. 
33 The University of Hong Kong, Pokfulam, Hong Kong 
34 Tsinghua University, Beijing 100084, P- R- China 
35 Universitaet Giessen, 35392 Giessen, Germany 
36 University of Hawaii, Honolulu, Hawaii 96822, USA 
37 University of Minnesota, Minneapolis, Minnesota 55455, USA 
38 University of Science and Technology of China, Hefei 230026, P. R. China 
39 University of the Punjab, Lahore- 54590, Pakistan 

40 University of Turin and INFN, Turin, Italy 
41 Wuhan University, Wuhan 430072, P. R. China 
42 Zhejiang University, Hangzhou 310027, P. R. China 
43 Zhengzhou University, Zhengzhou 450001, P. R. China 



b 



a also at the Moscow Institute of Physics and Technology, Moscow, Russia 
on leave from the Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 
c also at the PNPI, Gatchina, Russia 
d now at Nagoya University, Nagoya, Japan 
(Dated: December 2, 2011) 

The processes n' c — > p°p°, K*°K*°, and <j>4> are searched for using a sample of 1.06 x 10 8 ip' 
events collected with the BESIII detector at the BEPCII collider. No signals are observed in any 
of the three final states. The upper limits on the decay branching fractions are determined to be 
B(n' c -> p°p°) < 3.1 x 10" 3 , B{n' c -> K*°K*°) < 5.4 x 10" 3 , and B(n' c -> H) < 2.0 x 10" 3 at the 
90% confidence level. The upper limits are lower than the existing theoretical predictions. 

PACS numbers: 14.40.Pq, 12.38. Qk, 13.20.Gd, 13.25. Gv 



The radially (n=2) excited S-w&ve spin-singlet char- 
monium state, rj' c , labeled rj c (2S), was observed in 
B ± -> K±ri' c , r)' c -> K^K+tt* by the Belle Collabo- 
ration [l| and was confirmed by the CLEO and BaBar 
collaborations [2]. In addition to the KKir final state, 
rf c -> 3(vr+7r-), K + K-2(tt + tt-), KIK^-k+tt- , and 
7r + 7r~ K + K~n° are also reported [3|. The produc- 
tion of r/' c is also expected from the radiative magnetic 
dipole (Ml) transition of ip'. The decay ip' iv'o 
n' c -t KgK + Tr~ + c.c. was observed at BESIII [| with a 



branching fraction B(ip' -> ~/n' c ) = (4.7±0.9±3.0) x 10" 4 , 
confirming the possibility to study r/' c properties in ip' 
transitions. In this analysis, we search for the n' c decay- 
ing into vector meson pairs. 

The decay modes rf c — > VV, where V stands for a 
light vector meson, are supposed to be highly suppressed 
by the helicity selection rule @. But in Ref. a 
higher production rate of rj' c — > VV is predicted, tak- 
ing into consideration significant contributions from in- 
termediate charmed meson loops, which provide a mech- 
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anism to evade helicity selection rule Q • The intermedi- 
ate charmed meson loops can also significantly suppress 
if)' —> VP (where P stands for a pseudoscalar meson) 
strong decay amplitudes [§], which may help to explain 
the "/57T puzzle" in charmonium decays [9j . The measure- 
ment of B(r\' c — > VV) may help in understanding the role 
played by charmed meson loops in r\ c — > VV. 

In this study, an e + e~ annihilation data sample with 
(1.06 ± 0.04) x 10 s ip' events pjj is analyzed. Another 
data sample of 923 pb" 1 at \fs = 3.773 GeV is used to 
estimate non-ip' background. The data were collected 
with the BESIII detector which is described in detail 
elsewhere [ll| . A charged-particle tracking system, main 
drift chamber, is immersed in a 1 T magnetic field. A 
time-of-flight system and an electromagnetic calorime- 
ter (EMC) surrounding the tracking system are used to 
identify charged particles and to measure neutral particle 
energies, respectively. Located outside the EMC, a muon 
chamber is used to detect muon tracks. 

A Monte Carlo (MC) simulation is used to determine 
the mass resolution and detection efficiency, as well as to 
study backgrounds. The simulation of the BESIII detec- 
tor is based on GEANT4 [l2j, where the interactions of 
particles with the detector material are simulated. We 
use the program LUNDCRM [13| to generate inclusive MC 
events for the background study, where the branching 
fractions for known decay channels are taken from the 
Particle Data Group (PDG) (l4j |. For the signal channel 
ip' — > jt]' c , the photon is generated with the polar angle 
distribution 1 + cos 2 8. To generate the correct decay an- 
gle distributions, the rj' c — > VV decays are modeled with 
S VV model [l5| , and V decays are generated by the VSS 
model [HI, which is used to describe decays of a vector 
particle into two scalars. 

We search for the r)' c in three exclusive decay chan- 
nels: ip' -> 1P°P° -> 72(7r+7T-), V' jK*°K*° -> 
jtt + tt~K + K~, and ip' -> jcfxp -t j2(K + K~). These 
final states, denoted as ip' — > jX hereafter, contain one 
radiative photon and four charged tracks. The charged 
tracks are required to pass within 1 cm of the e + e~ an- 
nihilation interaction point transverse to the beam line 
and within 10 cm of the interaction point along the beam 
axis. Each track should have good quality in track fit- 
ting and satisfy |cos#| < 0.93, where 9 is the polar 
angle with respect to the e + beam direction. Recon- 
structed events are required to have four charged tracks 
and zero net charge. Information from dE/dx and time- 
of-flight is used for charged-particle identification (PID), 
and Xpid (*) i s calculated for each charged track, where 
i is the corresponding charged-particle hypothesis includ- 
ing pion, kaon, and proton. For a specific decay channel, 
the total Xpid ^ s obtained by summing Xpid W over 
the charged tracks. There is a loop to match the charged 
tracks to the final state particles in the decay channel, 
and the matching with the minimum Xpid 1S adopted. 
The decay channel for a reconstructed event is selected as 
the one with the minimum Xpid among possible decay 
channels. Photons are reconstructed by clustering EMC 



crystal energies with a minimum energy of 25 MeV. The 
photon candidates are required to be detected in the ac- 
tive area of the EMC (| cos# 7 | < 0.8 for the barrel and 
0.86 < | cos(9 7 | < 0.92 for the endcaps). Timing require- 
ments are used in the EMC to suppress electronic noise 
and energy deposits unrelated to the event. 

In order to reduce background from non-VV pro- 
duction, the invariant masses of the final decay parti- 
cles are required to satisfy 0.67 GeV/c 2 < AI w + n - < 
0.87 GeV/c 2 , 0.85 GeV/c 2 < M n ± K T < 0.95 GeV/c 2 , 
and 1.01 GeV/c 2 < M K+K - < 1.03 GeV/c 2 , for p°, K*° 
and 4> candidates, respectively, which are determined by 
fitting their mass distributions in the XcJ mass region. 
Here the background level has been considered in the 
choice of the selection criterion for each channel. The 
ratios of signal over non-V background are near 1 at the 
edges of the mass selection region for p° and K*°. 

A kinematic fit is performed to improve the mass reso- 
lution and reject backgrounds. The four-momenta of the 
charged tracks and the photon candidate are constrained 
to the initial ip' four- momentum (4C fit). When there is 
more than one photon, the photon with the minimum x 2 
from the 4C fit, xtc> ^ s taken as the radiative photon, 
and x\c * s required to be less than 40. 

Background from ip' — » TT + ir~J/ip with J/ip decay- 
ing into a lepton pair is removed by requiring the re- 
coil mass [l7| of any 7r + 7r~ pair to be below the J/tp 
mass (m^I < 3.05 GeV/c 2 ). Events from ip' -> r)J/ip, 
with t] — > 7r + 7T~7r° (7) and J /ip decays into lepton pairs, 
are also removed by this requirement. 

The background remaining can be separated into three 
categories: events with no radiative photon (ip' — > X); 
events with an extra photon in the final state (ip' — > ir°X, 
ti-O _ > ^); and events with the same final state as the 
signal (ip' — > jX), but where the photon comes from 
initial state radiation or final state radiation (FSR). 

The background from ip' — » X with no radiative pho- 
ton comes from events where the charged tracks plus 
a fake photon satisfy the 4C kinematic fit. In the 
X mass spectrum from a 4C kinematic fit, this back- 
ground contributes a peak close to the r)' c mass, around 
3.656 GeV/c 2 , and decreases sharply at high mass due 
to the 25 MeV requirement on the photon energy. If 
the measured energy of the candidate photon is not used 
in the kinematic fit, thus becoming a 3C fit, this back- 
ground lies around the ip' mass region (3.66 GeV/c 2 ~ 
3.70 GeV/c 2 ) in the mass spectrum, as the photon en- 
ergy from the fit tends to be close to zero energy (see 
Fig. Q}. There is little change in the r/' c mass resolution 
due to one less constraint in the kinematic fit, but the 
separation of the rj' c signal from the background is much 
improved. Therefore, the result from the 3C fit (Af| ■ ) is 
taken as the final mass spectrum. 

The background from ip' — > ir a X is measured from 
data by reconstructing the ir° from its decay into two 
photons. If there are more than two photons, the ir° 
candidate is selected as the one with the minimum x 2 
from a 5C fit (4C plus a ir° mass constraint). x\c < 30 
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. 3C for T| c (2S)->p°p° 
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FIG. 1: Comparison between 3C and 4C kinematic fits (un- 
normalized) . Shown in the plot are the signal with the 3C fit 
(filled circles), signal with the 4C fit (open circles), ip' — >■ X 
background with the 3C fit (solid line), and ip' X back- 
ground with the 4C fit (dashed line). 



is required to veto backgrounds. A MC sample of if)' — > 
tt°X is used to determine the efficiency ratio between 
events passing the ip' — > 7A and ip' — > ir°X selections. 
Finally, the efficiency ratio is used to scale the ip' — > 
n°X sample selected from data to obtain the background 
contamination from ip' — » tt°X as a function of the X 
invariant mass. This b ackg round, which is described with 
a Novosibirsk function [18j as shown in Fig. [2j contributes 
a smooth component in the Xcj (J = 0, 1,2) mass region 
(3.35 GeV/c 2 - 3.60 GeV/c 2 ), and is almost negligible 
above 3.60 GeV/c 2 . 

The background shape from ip' — > (jfsr)X is obtained 
from MC simulation, where the FSR photon is simulated 
with PHOTOS 0. The fraction of events with FSR is 

defined as Rfsr = ^ffj** , where N 7FSR x (Nx) is the 
number of events containing an (no) FSR photon that 
survive selection. This fraction is obtained from measur- 
ing the FSR contribution in ip' — ► 7Xc0> Xco — > {jfsr)X. 
The event selection of this FSR sample is very similar 
to that of the signal mode, except that the reconstructed 
final state contains two photons, where the softer pho- 
ton is regarded as the FSR photon. The energy of the 
FSR photon is not used when performing the 3C kine- 
matic fit for this sample. Events from ip' — > tt°X are the 
main background for the FSR sample and are excluded 
by requiring the invariant mass of the two photons to 
be outside of the 7r° signal region. Figure shows the 
two-dimensional distribution of M 5- c versus y . If 

we add the four-momenta of the FSR photon and X to 
calculate the invariant mass for events with M x below 
the X c0 mass in the PDG (Af£f G ), M^ srX peaks at 
Mjfg indicating the photon is indeed from FSR. As a 
result, events from Xco ~^ X are in the dashed-line box in 
Fig. [3J while events from Xco ~^ IfsrX are in the solid- 
line box in Fig. |3J In this way, we can obtain Rfsr for 



MC simulation and data. The factor Jfsr is defined as 
the ratio of Rfsr measured in data to that determined 
in MC simulation. This FSR measurement is performed 
for two final states; f F SR = 1-70 ± 0.10 and 1.39 ± 0.08 
are determined for X = 2(ir + Tr~) and X = tt + tt~ K + K~ , 
respectively. The errors are the statistical errors of the 
sample and the uncertainties of the background estima- 
tion. These factors are used to scale fractions of FSR 
background events [ip' —> ("fFSR)X} in the MC samples 
to estimate the background in data. 

Data taken at s/s = 3.773 GeV are used to esti- 
mate backgrounds from the continuum [e + e~ — > 7* — > 
{jfsr)X] and initial state radiation (e + e~ — > 'JisrX). 
MC simulation indicates that ip" decays contribute neg- 
ligible background in the modes under study. Using the 
luminosity normalization and energy dependence of the 
cross section, there are 46±3 and 8±2 background events 
expected for V = p° and V = K*°, respectively. For 
V = (p, no events survive the selection. 

The signal yields are extracted from an unbinned max- 
imum likelihood fit to the Myy distribution. The sig- 
nal shape is obtained from MC simulation, following 
_BVF(mo, r) x E^x damping , where mo and T are the mass 
and width of the Breit-Wigner for signal and XcJ, £7 is 
the cube of the radiative photon energy, which is nec- 
essary in an El/ Ml radiative transition, and damping 
stands for a damping function used to damp the diverging 
tail caused by the at lower mass region (correspond- 
ing to a higher energy radiative photon). One damping 

function used by KEDR [20( is defined as E Eq+ ^ -e ) 2 > 
where Eq is the most probable energy of the transition 
photon. It is also necessary to convolute this with a Gaus- 
sian function G(/x, a) to take the mass resolution differ- 
ence between MC simulation and data into account. The 
mean (//) and standard deviation (a) are free parameters 
for the XcJ signals. For r]' c , they are fixed to the values 
extrapolated from Xcj with a linear assumption. In the 
fit, the estimated backgrounds from ip' — > tt°X and the 
continuum are fixed. The shape of the ip' — > ("Jfsr)X 
background comes from the MC simulation. The fraction 
of MC data with an FSR photon is scaled by the factor 
fFSR to estimate the fraction of data with FSR back- 
ground. Figure 0] shows the final fitting results to the 
3C mass spectrum. The values of x 2 / n <if are 0.68 and 
0.72 for p°p° and K*°K*°, respectively, indicating good 
fits. The numbers of n' c events obtained are 6.5 ± 6.4 and 
6.9 ±4.8 for V — p° and K*°, respectively. No fit is per- 
formed for (pep, since there is only one rj' c — > (jxp candidate 
event in the signal region. 

The systematic uncertainties related to tracking, pho- 
ton reconstruction, PID and the kinematic fit are esti- 
mated with specially selected control samples [2l|. An 
efficiency can be defined as the ratio of Xcj yield for VV 
with the V mass requirement to that without this require- 
ment. The exact same method is applied to MC and the 
difference in the efficiency between MC simulation and 
data is taken as the corresponding systematic uncertainty 
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FIG. 3: The two-dimensional plots of Mjp versus x for events passing the ip' — > 'y^fshX selection with X — 2(it + tt~). 

From left to right they are (a) MC simulated XcO signal, (b) inclusive MC, and (c) data. In each plot the dashed-line and the 
solid-line boxes contain events without and with a FSR photon, respectively. MC simulations reproduce the shape well but 
not the amount of FSR events. 



caused by the V mass requirement, with the statistical er- 
ror included. An alternative damping function was used 
by CLEO (HI, exp(-£^/(8/3 2 )), which is inspired by the 
overlap of wave functions, with j3 = 65.0 ± 2.5 MeV from 
fitting the J jip — > "fr] c photon spectrum. The difference 
caused by the two damping functions is taken as a sys- 
tematic uncertainty. The main backgrounds that may 
affect our fit result in the rj' c mass region are the contri- 
butions from FSR in ip' — > jfsrX and from the con- 
tinuum. Therefore, the systematic uncertainty from the 
background shape is estimated by changing the FSR and 
continuum contributions by la. There are also system- 
atic uncertainties related to the mass and width of the T]' c , 
which are estimated by comparing the r)' c yields with the 
mass and width fixed to the center values or randomly 
selected values according to a Gaussian distribution. Ta- 
ble[T]shows a summary of all the systematic uncertainties. 

As there is no significant i]' c signal in any of the three fi- 
nal states, we determine upper limits on the ip' — ¥ 777^, — ¥ 
jVV production rates. We assume all the signal events 
from the fit are due to rj' c — > VV, neglecting possible 



TABLE I: The systematic uncertainties in the measured prod- 
uct branching fraction B(ip' — > "/rj' c ) x B(rj' c — > VV). 



Source 


P° 


K*° 




Background (%) 


14.9 


9.9 


0.0 


Tracking (%) 


8.0 


8.0 


8.0 


Photon reconstruction (%) 


1.0 


1.0 


1.0 


Particle ID (%) 


8.0 


8.0 


8.0 


4C fit (x 2 selection) (%) 


4.0 


4.0 


4.0 


V mass selection requirement (%) 


2.6 


1.1 


1.6 


Damping function (%) 


40.5 


10.0 


0.0 


Mass and width of rj' c (%) 


6.6 


5.8 


0.0 


Number of i>' (%) 


4.0 


4.0 


4.0 


Total (%) 


45.6 


19.9 


12.8 



interference between the signal and nonresonant contri- 
butions. The probability density function (PDF) for the 
expected number of signal events is smeared with the sys- 
tematic uncertainties (by convolution). For V = p° and 
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10" E| 




-1 I . ,-,-i-T i i i I i i i i I . I . . . . I /. . . Iili . u . II 

3.35 3.4 3.45 3.5 3.55 3.6 3.65 3.7 

M^ p „(GeV/c 2 ) 




FIG. 4: Invariant mass distributions of the vector meson pairs 
after a 3C kinematic fit for the modes (a) p°p°, (b) K*°K*°, 
and (c) (f>tp. Dots with error bars are data, and the solid curves 
in (a) and (b) are from the best fit to the mass spectra. No 
fit is performed for (c) due to low statistics. In (a) and (b), 
the T]' c signals are shown as short dashed lines, ip' — > 7r X 
backgrounds are in dotted lines, continuum in long dashed 
lines, and ip' — > (■Jfsr)X in short dash-dot-dotted lines. 



K*°, the PDF is taken to be the likelihood distribution 
in fitting the invariant mass distributions in Fig. |4] by 
setting the number of rf c signal events from zero up to 
a very large number. For V = <f>, the one event in the 
r)' c mass region is taken as signal for simplicity, and the 
PDF is assumed to be a Poisson distribution. 



The upper limit on the number of events 
at the 90% C.L., N* vv , corresponds to 

N up 

J yVV PDF(x)dx/ / °° ¥T)¥{x)dx = 0.90 on the smeared 
PDF. The left half of Table [CD shows N u p, the efficiencies 
from MC simulation, and the upper limits on the 
product branching fraction B(i/j' — > jri' c ) x B(r]' c — > VV). 
Using B(tf>' -> 777^) = (4.7 ± 0.9 ± 3.0) x 10" 4 Q, the 
corresponding upper limits on B(r/' c — > VV) are listed in 
the right half of Table HI In calculating B up (t]' c -> VV), 
the error on B(ip' —> 777^.) i g taken as a systematic 
uncertainty to smear the PDF. The theoretical predic- 
tions on branching fractions for r]' c — > VV, which are 
calculated with T n > = 10.4 ±4.2 MeV HH, are also listed 
in Table HU 



In conclusion, no obvious r)' c signal was observed in 
decays into vector meson pairs: p°p°, K*°K*°, and 
<t>(j). The upper limits on the product branching frac- 
tion B(ip' — > jri' c ) x B(r/' c — > VV) and r]' c decay branching 
fraction B(r)' c — > VV) are determined. These upper lim- 
its are smaller than the lower bounds of the theoretical 
predictions @, although the difference is very small for 
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TABLE II: From left to right, they are efficiency, upper limits at the 90% C.L. on the yield, product branching fraction 
B(ip' — > 77/c) x 13{v'c VV), rj' c decay branching fraction B(r/' c -> VV), and theoretical predictions from Ref. [g]. 



V 


e (%) N^ v - 


-> Wc -► ivv) (10- 7 ) 




-> VV) (10 -3 ) 




p° 


14.3 19.2 


12.7 




3.1 


6.4 to 28.9 




16.5 15.2 


19.6 




5.4 


7.9 to 35.8 


4> 


19.9 3.9 


7.8 




2.0 


2.1 to 9.8 
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